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Abstract 

A personal computer-based system was recently prototyped that uses measurements from a three-axis 
magnetometer (TAM) to estimate the attitude and rates of a spacecraft using no a priori knowledge of the spacecraft’s 
gats. P asIsM ie; using in-flight data from the Sola,. Anomalous, and MagneiostLoc Pafficl^Z- S? oh 

i «; ° thS , S ^® m , and demons,ra,ed ,hat attitude and rate estimates could be obtained accurately to 

1.5 degrees (deg) and 0.01 deg per second (deg/sec), respectively, despite limitations in the data and in the 

Sn C fTn fnlh ? p rth tP? P o Per StUd ' eS the accurac V of the Kalma " fitter in the system using several 

orbits of in-flight Earth Radiation Budget Satellite (ERBS) data and attitude and rate truth models obtained from hiqh- 
precision sensors to demonstrate the practical capabilities. ” 

This paper shows the following: 

* telemetered TAM data, attitude accuracies of 0.2 to 0.4 deg and rate accuracies of 0 002 to 
minlmah^n^g of'the fitter^ atMUde contro1 rec l uirem ents of 1 deg and 0.005 deg/sec) can be obtained with 

• ^Placing the TAM data in the telemetry with simulated TAM data yields corresponding accuracies of 0 1 to 

0.2 deg and 0.002 to 0.005 deg/sec, thus demonstrating that the filter's accuracy can be significantly 
enhanced by further calibrating the TAM. y y 

Factors affecting the filter’s accuracy and techniques for tuning the system's Kalman filter are also presented in this 
paper. 


1. Introduction 

Contingency algorithms for determining the attitude and rates of a spacecraft using only a three-axis magnetometer 
(TAM) and using no a priori knowledge of the spacecraft’s state have been of interest recently in view of situations 
such as the control anomaly on the Earth Radiation Budget Satellite (ERBS) in 1987 (Reference 1) and the control 
system failure on the Relay Mirror Experiment (RME) satellite (Reference 2). It has recently been shown (Reference 
3) using in-flight data from the Solar, Anomalous, and Magnetospheric Particles Explorer (SAMPEX) mission that 
an optimal combination of two different algorithms — the Deterministic Attitude Determination From Magnetometer- 
Only Data (DADMOD) (Reference 2) and the Real-Time Sequential Filter (RTSF) (References 4 and 5)— solves the 
problem by maintaining both accuracy and robustness. The concept was subsequently automated by prototyping a 
personal computer (PC)-based system (Reference 6) referred to here as TAMONLY for the Flight Dynamics Division 
(FDD) at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC). 

The working principle of TAMONLY involves using DADMOD to obtain coarse estimates of the epoch state of the 
spacecraft and then using this state as the initial state for the RTSF’s finer estimates. The previous study of 
TAMONLY (Reference 6) involved extensive tests using a variety of SAMPEX in-flight data and demonstrated that 
accuracies of 1 .5 degrees (deg) in attitude and 0.01 deg per second (deg/sec) in spacecraft rates could be obtained. 
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Note that these accuracies are within the attitude requirements for SAMPEX (2 deg). These accuracy estimates, 
however, are only coarse because of the following factors: 

• The most accurate sensor onboard SAMPEX is the Sun sensor, which is accurate to 0.5 deg. 

• SAMPEX does not have gyros so that the truth model for the rates was obtained by differentiating the attitude 
truth model data. 

• Short spans of data (< 3000 sec) were used. 

• The TAM was not calibrated. 

• The TAM is used to determine the pitch angle of the attitude truth model. 

TAMONLY was also used to study the above-mentioned ERBS control anomaly, and similar accuracies were 
obtained (Reference 7). Here again, accurate truth models were absent and only coarse accuracy estimates could be 
obtained during periods when gyro data were unsaturated in the B-Dot control mode of the spacecraft. 

This paper attempts to remedy the situation by focusing on accuracy evaluation of the RTSF in TAMONLY using 
ERBS in-flight data The data are from a period when the spacecraft was in its nominal mission mode soon after its 
launch in 1984; during this period, accurate truth models are available for the attitude via Earth-sensor and fine- 
Sun-sensor measurements, and accurate gyro data provide the truth model for the rates. 

Section 2 of this paper describes the Kalman filter, and Section 3 describes ERBS and the characteristics of the data 
used in this study. A discussion of the criteria used to evaluate the RTSF is given in Section 4. The results obtained 
using telemetered TAM data and the results obtained using simulated TAM data are given in Sections 5 and , 
respectively. Section 7 presents a summary of the conclusions of the study. 

2. Description of the Kalman Filter 

In view of space considerations, only details relevant to the tuning of the RTSF are presented here. A full 
mathematical description of the RTSF has been provided elsewhere (References 4 and 5). 

The RTSF’s state vector X comprises of the four components of the attitude quaternion, q , and the corrections, b , 
to the spacecraft’s rates, © : 

X = [q T b T f 0) 

(Note that the components of b and to are resolved along the spacecraft s x, y, z axes.) 

The RTSF uses sensor data to estimate q as well as b , with b being estimated kinematically in the same manner 
as gyro biases for a gyro-based spacecraft; i.e., by attributing differences between the measured and propagated 
attitudes to errors in <n. The b estimates are then used to correct d>, and these corrected rates are used as initial 
conditions to propagate Euler’s equation to the next measurement time. The propagation of b is modeled via a first- 
order Markov model: 


where is a white noise term, and x is a finite time constant. A suitable value for x is the time between 
measurements -5 seconds for the SAMPEX data and 16 seconds for the ERBS data used here. (In contrast, the same 
model, when used for gyro bias estimation, requires ax of several hours.) 

The rates are assumed to contain a white noise component, f) a , and are propagated using Euler s equation after 
accounting for the angular momentum contributed by the wheels, and the total external torques acting on the 
spacecraft. TAMONLY currently models the gravity-gradient torque and the magnetic control torque acting on the 
spacecraft. (The aerodynamic drag torque and the radiation pressure torque have been intentionally omitted to 
reduce the amount of spacecraft modeling required. The RTSF relies on the rate-corrections, b , to compensate for 
the small effects of these torques.) 
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The covariance matrix, P, is propagated by numerically integrating the following equation. 

~ = F(^)P+PF T (a,) + Q (3) 

Here F(a>) is described in Reference 5; the quantity of interest is the 6 x 6 matrix Q that quantifies the propagation 
noise and is of the following diagonal form: 


Q = diag [ Q ax , Q uv , Q av Q hx , Q hyt Q hz ] 


(4) 


Here Q ax , Q ay , and Q ai represent the strength of the noise term rj„ and contribute to the growth of the attitude error 
covariances about the body X-, Y-, and Z-axes during propagation. Similarly, Q bx , Q by , and Q bz represent the 
s rength of the noise term r \ b , and contribute to the growth of the error covariances of the rate corrections during 
propagation. These six elements of Q are collectively referred to here as the Q-parameters and rough estimates of 
t eir numerical values for ERBS are given in the next section. Another quantity that we must consider during 
tuning is o TA m, the strength of the white noise in the TAM measurements. 


3. Description of ERBS and Data Characteristics 

ERBS is an Earth (geodetic)-pointing three-axis stabilized spacecraft with nominal inclination 57 deg and a nominal 
altitude of 610 kilometers (References 8 and 9). The onboard attitude requirement of 1.0 deg is met by using a fine 
Sun sensor of accuracy 0.05 deg and a scan wheel-type Earth sensor of accuracy 0.02 deg. These high-precision 
sensors are complemented by a TAM onboard the spacecraft. Rate information is obtained from gyros that are 
accurate to 0.001 deg/sec. Attitude control is achieved through magnetic torquer bars (MTBs) and a momentum 
wheel that provides momentum bias along the pitch axis (negative orbit normal). A small angular momentum 
contribution also arises from the scan wheels in the Earth sensor assembly. ERBS also uses thrusters for its monthly 
180-deg yaw maneuvers, but these were not activated during the data period discussed here. 

The data used in this study are a 13-hour span of processed (adjusted) engineering data generated on the mainframe 
computers using the ERBS Attitude Ground Support System (AGSS) and covering the period 841225.221932 to 
841226.105348. The adjustments comprised mainly bias determinations for the Sun and Earth sensors, the gyros, 
and the TAM. It is important to note here that the TAM was not calibrated for misalignments and MTB coupling 
As shown in Reference 10, these calibrations can significantly alter the performance. The attitude history generated 
by the AGSS using the adjusted Earth and Sun sensor data is taken as the attitude truth model, and the adjusted gyro 
data are taken as the truth model for the rates. For the remainder of this paper, the times of events are described 
relative to the epoch of 841225.221932. 

These data were then downloaded to a PC where a dataset having the TAMONLY input format was generated. The 
processed telemetry was input to TAMONLY at a period of 16 sec (equal to the availability of control torque data) 
and consisted of the following: 

• Time 

• Reference magnetic field 

• Measured magnetic field 

• Dipole moments of the MTBs 

• Wheel speeds 

• Position and velocity of the spacecraft 

The profile of the truth models is presented in Figures la through lc, where Figure la shows the true attitude. Figure 
lb shows the true rates (where the rate components in the figure are denoted by “wx” “wy” “wz”), and Figure lc 
shows the TAM residuals generated by the true attitude. Note the statistics of the data presented in Figures la 
through lc. For example, it can be seen that the mean residuals are different from zero; in fact, these mean values 
had to be subtracted off from the TAM measurements input to TAMONLY to improve its performance slightly 
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Figure la. Profile of the Attitude Truth Model Obtained Using Earth and Sun Sensors 
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Figure 1b. Profile of the Rate Truth Model Obtained Using Gyros 
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Figure 1c. TAM Residuals Obtained Using the Attitude Truth Model 


During the course of this investigation, the optimal tuning parameters for the telemetered data (Section 5) indicated 
that calibration of the TAM might be in order. To see what effects a perfect TAM might have, simulated TAM data 
with no noise were generated using the true attitude and the reference magnetic field. These simulated TAM data 
were then used m place of the actual TAM measurements thus providing a much-needed degree of control over the 

rf S flf descnbe ? Section 6 Note that only the TAM data were replaced; other quantities, such 
as the MTB and wheel data, were left untouched in this process. 

The noise parameters were determined as follows by considering the noise in the data. The nominal value of 
was determined as 6.4 milligauss (mG), which is the discretization of the TAM data in the telemetry. Regarding the 
^'Parameters it was found that the principal source of noise during propagation was the digitization in the 
telemetered wheel momenta, which was about 0.7 Newton-meter-second (N-m-s) along the Y-axis and 0 01 N-m-s 
Z ' ax,s ‘ In view of the teIemetr y P eriod of 16 sec, this generated torques of the order of 0.04 N-m and 
0 . 00 ". N '™’ respe 5 lvely - 11 was assu med that the noise torque along the X-axis was of the same magnitude as that 
about the Z-axis. The squares of the corresponding uncertainties in the spacecraft rates are then 

(Aco^) 2 = 2.3 x 10‘ 8 radians 2 /second 2 (rad 2 /sec 2 ) c<\ 


(Aco,) 2 = (Aco 2 ) 2 =10 " rad 2 /sec 2 (6) 

During the telemetry period. At, the attitude error covariances corresponding to the i ,h axis grow by (Aco ) 2 (A/1 2 rad 2 
Because, from Equation (3), Q is the rate of change of P, this yields: 

Qai = (Act),) 2 At rad 2 /sec, (i = x, y, z ) ( 7 ) 

where At = 16 sec and Equations (5) and (6) are to be used for (Aco,) 2 . 
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The bottom three elements of Q are assigned numerical values by invoking the following relationship (Reference 1 1 ) 
between the steady-state covariance, p„ of an exponentially correlated random variable (such as b ) and its spectral 

density, q: 

a = IEzl ( 8 ) 

T 

It is assumed that the noise in b is of the same order of magnitude as that in the rates themselves, so that p- was 
chosen equal to the squares of the rate uncertainties to yield 

Q hi = 2 (Ato,) 2 / x rad 2 /sec 3 (i = x, y, z) 

where x = 16 sec and the numerical values of (Ato, ) 2 are to be used from Equations (5) and (6). 

From Equations (7) and (9), it can be seen that the Q-parameters can be parametrized by the : rate : uncertainties 
according to this scheme, and these uncertainties will be used in describing the results of tuning the RTSF. It should 
be noted that this study is limited to tuning the RTSF noise parameters to improve its performance, and that the 
errors contributed by imperfect models of the spacecraft’s dynamics are ignored here. The latter contributions may, 
in fact, be sizable because a great number of error sources exist, such as spacecraft mass properties (inertia tensor), 
flexible modes of the spacecraft, misalignments of the wheel axes, scale factors for the wheel speeds, misalignments 
of the MTBs, and limitations of the numerical integration algorithms. 

4. Accuracy Criteria 

The following criteria were adopted to determine the accuracy of the RTSF as well as the optimality of the tuning 
parameters: 

• Attitude differences , which are the differences between the true attitude and RTSF estimates 

• Rate differences, which are the differences between the true rates and RTSF estimates 

• TAM residuals, which depend upon the accuracy of the propagation, the accuracy of the reference magnetic 
field, and the accuracy of the TAM. Accurate RTSF estimates are typified by TAM residuals of less than 

10 mG. 

• The RTSF’s rate corrections, b , which are not to be confused with the rate differences noted above. 
Accurate RTSF estimates are typified by rate-correction estimates of less than 2 deg/hour. In this study, the 
rate-correction estimates are used only qualitatively, i.e., they are used to confirm the convergence of the 
RTSF to the truth model. 

DADMOD usually initializes the RTSF to within 5 deg of the true attitude and 0.05 deg/sec of the true rates. The 
results from the first 5000 sec of a run were omitted from the statistics to allow for filter convergence. Typical 
results of a TAMONLY run are presented in Figure 2a and 2b where results for the Y- and Z-axes are shown. In 
these figures, the rate components are denoted by “wy” and “wz,” and the components of the RTSF’s rate corrections 

are denoted by “by” and “bz.” 

5. Results Using In-Flight ERBS Data 

About 25 runs were made using a wide range of the noise parameters. A notable (and welcome) feature was that a 
wide range of numerical values did not alter the accuracies substantially. A sample of the results is shown in 
Table 1. The first column of Table 1 denotes the numerical values of the noise parameters [a ww , (Acu y ) , and 
(Ato) 2 ], together with a mnemonical description of the set. Against each set of tuning parameters, the mean and 
root-mean-square (r.m.s.) of the attitude differences, rate differences, and TAM residuals for each body axis are 
displayed. Note that accuracies of the RTSF are taken to be the r.m.s. values of the differences from the truth 

model. 


184 


Error (deg/sec) wy Error (deg/sec) 


mean = -0.000002 degteec, r. m. s. = 0.0055 deg'sec 



Figure 2a. Typical Attitude Results (Top Two Plots) and Rate Results 
(Bottom Two Plots) Obtained Using TAMONLY 


185 






Z residual (mG) Y residual (mG) 



Time since 841225.221932 (sec) x 10* 


Figure 2b. TAM Residuals (Top Two Plots) and Rate-Correction Estimates 
(Bottom Two Plots) Obtained During the Run of Figure 2a 
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A number of useful deductions can be made by examining Table 1 : 


' <o - 2 ,o 04 *■ in - «« » 

• s ™ * ,he — rf 

* JSLSfas^js! TAM residuais comparab,c in - *- - «*- .* « 

* 2 Shc*" K f Jl" "uwTms TTI i"t atts a n “ d for calibra,ins ““ TAM fOT 

P~ than ° U, ” eri ‘ al ^ f ° f ** Q 

numte — ■ «*- * * 
(5864 sec) and (2) spanning 4000 sec each. Sarnie ^ T “* 

penod effects were detected in the statistics. P T b 2 d 2b ’ n ° s,gnificant orbital 

6. Results Using Simulated TAM Data 

l7^"sS?n "^“.,'trL^ replaced " i ‘ h si ™'^ TAM data and the 
parameter results of Table 1 with those of Table 3: UC '° nS ^ ** made by COmpanng the °P timal 

• resuiK - — - *» - — 

• The TAM residuals are significantly smaller than in Table 1 . correctly reflecting the perfect sensor situation. 

reSUllS <third row of Table 31 « obtained when the modeled noise parameters for ERBS are used 
and result accurac.es of 0. 1 0.2 deg in attitude. Note that, although 0,1 = 0 Sl7a mal 

,o avoid di — ■ n ° ----- 

• z - ~ - — 
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Table 2a. Error Statistics Over Orbital Period 
(Using Optimal Parameters of Table 1) 


Time Span 
(seconds) 

Attitude Accuracies (deg) 

Roll 

(mean, r.m.s) 

Pitch 

(mean, r.m.s) 

Yaw 

(mean, r.m.s) 

5008-10860 

-0.047, 0.150 

0.098. 0.470 

-0.092, 0.350 

16752-22640 

-0.025, 0.207 

-0.070, 0.345 

-0.106, 0.424 

34384-40256 

-0.130, 0.223 

-0.004, 0.233 

-0.357, 0.631 


Table 2b. Error Statistics Over 4000-Second Spans 
(Using Optimal Parameters of Table 1) 


Time Span 
(seconds) 

Attitude Accuracies (deg) 

Roll 

(mean, r.m.s) 

Pitch 

(mean, r.m.s) 

Yaw 

(mean, r.m.s) 

9008-13008 

-0.044, 0.169 

-0.025, 0.391 

-0.165, 0.338 

21008-25008 

-0.043, 0.200 

-0.171,0.187 

-0.329, 0.480 

37008-41008 

-0.213, 0.296 

-0.008, 0.274 

-0.282, 0.373 


Table 3. Error Statistics and Tuning Using Simulated TAM Data 


Description 
[OTAMh (Aav) 2 * (A(Dj) ] 

Attitude Difference* With Truth 
(deg) 

Mean / r.m.». 

Rate Differences With Truth 
(deg/tec) 

Mean / r.m.s. 

TAM 

Residuals 
Mean / r.m.i 

(mG) 

t. 

x 

y 

Z 

X 

y 

Z 

X 

y 

z 

Optimal Parameters 
of Table 1 

r a A 3*1 O'® 3x1 O' 8 ! 

-0.134/ 

0.207 

0.042/ 

0.313 

-0.260 / 
0.373 

0.0026 / 
0.0029 

0.00015/ 

0.0055 

0.00060/ 

0.0016 

-0.312/ 

1.055 

0.053 / 
0.574 

-0.301 / 
0.774 

l4Uf IU 1 wAIV J 

Small o tam 

/4A-2 Qvin"® 3*10'®! 

-0.117/ 

0.230 

0.003/ 

0.100 

-0.190/ 

0.282 

0.0026 / 
0.0031 

0.00010/ 

0.0054 

0.00089 / 
0.0019 


0.062/ 

0.269 

-0.306 / 
0.423 

(lU , JXIU , ja i u ; 

Modeled Parameters 
(10~ z , 2.3x10 ®, 10’ 11 ) 

-0.093/ 

0.179 

0.0004/ 

0.088 

-0.137/ 

0.219 

0.0026 / 
0.0030 

0.00008/ 

0.0053 

0.001 / 
0.0019 


0.046/ 

0.255 

-0.291 / 
0.414 


7. Conclusions 

The principal conclusions derived from this study are as follows: 

. Using telemetered TAM data, the accuracies (r.tn.s. diffemnces from the truth models) of TAMONLY were 
found to be 


0.2 to 0.4 deg for attitude 
0 002 to 0.005 deg/sec for rates 


rates. 
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• The results are not severely degraded if the default (SAMPEX) parameters are used. This is a very useful 
result, since tuning a Kalman filter is a manpower-intensive exercise. 

• No significant orbital-period effects were seen in the error statistics. 

• The simulated ERBS TAM data were useful by demonstrating the following: 

Additional calibration of the TAM for misalignments and MTB coupling effects will extract better 
performance from the RTSF. 

- The modeled noise parameters indeed give better results when the TAM data are of high quality. 

- Very little improvement resulted in the r.m.s. values for rate differences, thus indicating that further 
improvements in accuracy must be generated by more accurate dynamics, i.e., through improved 
propagation algorithms and more accurate spacecraft and torque models. 
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